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ABSTRACT 

We use the combined data-sets of the Millennium I and II cosmological simulations to 
revisit the impact of mergers in the growth of bulges in central galaxies in the ACDM 
scenario. We seed galaxies within the growing CDM haloes using semi-empirical rela- 
tions to assign stellar and gaseous masses, and an analytic treatment to estimate the 
transfer of stellar mass to the bulge of the remnant after a galaxy merger. We find 
that this model roughly reproduces the observed correlation between the bulge-to- 
total (B/T) mass ratio and stellar mass (M*) in present-day central galaxies as well 
as their observed demographics, although low-mass B/T < 0.1 (bulgeless) galaxies 
might be scarce relative to the observed abundance. In our merger-driven scenario, 
bulges have a composite stellar population made of (i) stars acquired from infalling 
satellites, (ii) stars transferred from the primary disc due to merger-induced perturba- 
tions, and (iii) newly formed stars in starbursts triggered by mergers. We find that the 
first two are the main channels of mass assembly, with the first one being dominant for 
massive galaxies, creating large bulges with different stellar populations than those of 
the inner discs, while the second is dominant for intermediate/low-mass galaxies and 
creates small bulges with similar stellar populations to the inner discs. We associate 
the dominion of the first (second) channel to classical (pseudo) bulges, and compare 
the predicted fractions to observations. We emphasize that our treatment does not 
include other mechanisms of bulge growth such as intrinsic secular processes in the 
disc or misaligned gas accretion. Interestingly, we find that the evolution of the stellar 
and gaseous contents of the satellite as it spirals towards the central galaxy is a key 
ingredient in setting the morphology of the remnant galaxy, and that a good match to 
the observed bulge demographics occurs when this evolution proceeds closely to that 
of the central galaxy. 

Key words: galaxies: formation galaxies: evolution galaxies: bulges galaxies: 
interactions galaxies: structure. 



1 INTRODUCTION 

In the A Cold Dark Matter (ACDM) cosmogony, galactic 
discs generically form inside the growing CDM haloes, while 
spheroids (bulges and elliptical galaxies) are thought to grow 
through mergers, by disc secular intern al processes, a nd/or 
by mis-aligned/perturbed infalling gas (|Mo et alj|20ld . and 
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references therein). Several arguments collected over the 
past decades suggest that major mergers of galactic discs 
produce prominent spheroids that are gravitationally sup- 
ported by random motions. Spheroids built this way might 



bulge-dominated galaxies (e.|. 
Gerhard! Il98l|; ISchweizerl 11982 



Toomre 1 19771: White 


1978 


Neeroponte & White 


1983 



Barnes! 19881 : lHerncmistlll992l ). The demographics of these 
(early-type) galaxies are thus expected to be tightly con- 
nected to their past galaxy merger activity, which is linked to 
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the major merger history of their haloes. An open question is 
whether the observed morphological distribution of galaxies, 
globally measured by the bulge-to-total mass (B/T) ratio, 
and its evolution with redshift, is consistent with the halo 
merger rates predicted by the ACDM cosmogony. 

In order to answer this question, it is necessary to un- 
derstand how tight is the relation between the growth of 
spheroids and the merger histories of their host haloes. The 
merger-driven growth of spheroids along the evolutionary 
path of galaxies is more complex than assuming that this 
growth is proportional to the history of halo mergers with a 
mass ratio above a given threshold. This caused by both, the 
non-linear relation between halo mass and s t ellar mass(e.g., 



Conrov fc Wechslerll20o"9l;lMoster et al.l2010l;lBehroozi et all 



2010; iFirmani fc Avila- Reese! 
Yang et al.l l201lJ ) 



20TB; 



Leauthaud et al] 



20121 : 



and the connection between 



progen- 
m ergers 



itor gas fractions and sphe r oid formation 
fe.g. iBarnes fc Hernquisdll996l: ISpringel fc Hernauistll2005l; 



et al] 



Governato et al. II2009I: iHammer et al.l l2009l : iHopki: 
2009bl: IStewart et al]|2009l ') . 

During a merger, the bulge in the primary galaxy may 
grow not only by acquiring stars from the merging satellite 
(the secondary), but also by violently relaxed stars trans- 
ferred from the primary disc, and by n e wly cre ated stars 
formed in starbursts. In [Hopkins et al] (|2009al . hereafter 
H09a), the authors developed analytic prescriptions to de- 
scribe these physical processes with the aim of predicting 
the amount of stellar material that is finally deposited into 
the bulge of the remnant. The parameters of these prescrip- 
tions were calibrated to roughly agree with the outcome of 
a large suite of full hydrodynamical simulations. 

Numerical s i mulat io ns (e.g . lAthanassoulal 120051 : 
lAvila-Reese et al] 120051 : ICombesi 12009] , and references 
therein) show that bulges can also grow through in- 
trinsic secular transport of angular momentum and 
dynamical heating of the stellar disc, giving rise to the 
so-called pseudobulges. In contrast with classical bulges, 
pseudobulges have more rotational support and share 
some of the properties of the inner discs (for rev iews see 
iKormendv fc Kennicutll2004 lusher fc DrorvteOOST ). Demo- 
graphics of bulge types in the very local environment show 
that low- mass galaxies have small B/T ratios (< 0.2) with 
spheroids that mostly belong to the pseudobulges category. 
This has been raised as a potential issue for ACDM which 
would seemingly predict higher fractions of merger-driven 
bulges, that ty pically have higher B/T ratios and are of the 
classical type (|Kormendv et alj|2010l; iFisher fc Drorvll20l"ll . 
hereafter FDll^ see also IWeinzirl et al] I2009T hereafter 
W+09). However, when the stellar mergers have a low mass 
ratio < 0.1 (see Section 3.2), the bulge may actually be 
predominantly populated by stars in the primary disc that 
are transferred to the centre and dynamically heated by 
instabilities induced by the merger. This is the dominant 
effect in this kind of merger rather than the addition of 
stars coming from the secondary. It is thus possible that 
bulges formed in minor/minuscule stellar mergers look like 
pseudobulges. 

Several studies have aimed to establish the connection 
between the halo merger history with the final galaxy B/T 
ratio (morphology). Some of these works are based on the 
semi-empirical halo occupation fr amework l|Stewart et al] 
120091 : iHopkins et alj|2009bl , l2010aT ). while others are based 



on semi-analytic models (SAMs) llKhochfar fc Silkl 120061: 
Parry et al. 2009; Benson fc Devereuxll2010l ; iDe Lucia et al] 



20 111 ; iFontanot et al] l201ll ). The most general results of 



these works are that: (i) the mapping of halo-halo mergers to 
stellar galaxy-galaxy mergers is far from linear and strongly 
depends on mass and redshift, (ii) the inclusion of the galaxy 
gas content in mergers significantly reduces the final B/T 
ratio, specially for low-mass galaxies and at higher redshifts, 
and (iii) the B/T ratio predicted in the ACDM scenario in- 
creases with stellar mass in a similar way as observations, 
although there seems to be fewer predicted bulgeless galaxies 
than observed. 

Our goal is to revisit the merger-driven bulge forma- 
tion in the context of the ACDM cosmogony and to com- 
pare the results with the observed bulge demographics. 
Here, we focus only on present-day central galaxies, al- 
though the population of satellites is considered along the 
evolution of central galaxies. The backbone of our model 
is an approach where, along the mass aggregation histories 
(MA Hs) of haloes taken f rom the Millennium Simula tions 
(MS, ISpringel et al.l 120051 : iBovlan-Kolchin et al.l 120091 ') . the 
galaxy stellar and gas masses are assigned according to to 
semi-empirically inferred average M*(Mh, z) and M 9 (M*, z) 
relations. The galaxy evolutionary tracks are quite different 
to the halo MAHs and are such that at each z there is a tran- 
sition stellar mass, Mt Ta ,n{z = 0) ~ 2x 10 10 Mq, above which 
the M» growth has halted and below which M» is actively 
growing (|Firmani fc Avil a- Reese 20T0I'). The former happens 
at earlier epochs for more massive galaxies (archaeological 
downsizing) and for the latter, less massive galaxies are more 
active at late epochs (downsizing in specific star formation 
rate). Note that our approach is different to that of SAMs 
since instead of modeling the complex galactic physics, the 
main galaxy properties are assigned at each epoch according 
to empirical relations. 

The paper is organised as follows. In Section [2] we de- 
scribe the way we estimate the central halo-halo merging 
time, and our semi-empirical approach to seed baryonic cen- 
tral galaxies into the growing progenitors of the halo popu- 
lation. In 5 ^3.11 we present the stellar merger fractions as a 
function of mass and redshift and compare them with recent 
observational estimates. In § £13.21 we give our main predic- 
tions related to: (i) the growth of bulges as a function of 
mass, and the contribution to this growth from the different 
mechanisms of bulge assembly; and (ii) the B/T demograph- 
ics of central galaxies at z = 0. In 5 ^3.2.31 our predictions 
are compared with current observational results. Finally, a 
summary and our conclusions are given in Section [4] 



2 SEMI-EMPIRICAL MODEL OF BULGE 
GROWTH 

2.1 N-body simulations 

We use the combined data-sets of the Millennium (MS-I) 
and Millennium II (MS-II) s imulations l|Springel et al.ll2005l ; 
IBovlan-Kolchin et al] [2009) that share the same particle 
number (2160 3 ) and were done on the context of a WMAP1 
cosmology with parameters: fl m = 0.25, 57a = 0.75, h — 
0.73, as = 0.9 and n s = 1; where fi m and Qa are the 
contribution from matter and cosmological constant to the 
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mass/energy density of the Universe, respectively, h is the 
dimensionless Hubble constant parameter at redshift zero, 
n s is the spectral index of the primordial power spectrum, 
and erg is the rms amplitude of linear mass fluctuations in 
8 ft -1 Mpc spheres at redshift zero. The MS-II has a box size 
(L — 137 Mpc on a side) that is 5 times smaller than the 
one of the MST, and thus it has a mass resolution limit 125 
times smaller: 9.4 x 10 6 Mq. Combining both simulations we 
can follow up to high redshift (z ~ 10) the merger and accre- 
tion histories of haloes having a wide mass range at z — 0: 
10 10 - 1O 15 M . 

We should note that the WMAP1 cosmological param- 
eters used in the simulations are different to those currently 
preferred by the 7-year WMAP results. In particular, as is 
lower and fi m is higher in WMAP7 which produces com- 
pensating effe cts in the abund ance and clustering of dark 
matter haloes. IGuo et ail (|2012T ) studied the impact of both 
cosmologies on galaxy formation using a SAM and found 
that both predict a similar evolution in the galaxy prop- 
erties since z — 3 with only a slightly lower autocorre- 
lation function at separations < 1 Mpc, particularly for 
lower masses, and also a lower fraction of satellites. Hence, 
a slightly lower merger rate between satellites and centrals 
(producing remnants with lower B/T ratios) would be ex- 
pected if a WMAP7 cosmology was used instead of the one 
we adopt here. The difference however is not expected to be 
significant. 



(i) the start of the merger (z s tart)> i- e -i when the sub- 
halo was part of an independent FOF halo for the last 
time. The halo mass ratio of the merger is defined by us- 
ing the subhalo and main halo masses at this time, jUhalo = 

Zstart J , 

(ii) the "dissolution" of the subhalo (tdiss), i.e., when the 
merged subhalo at time U can no longer be resolved as an 
independent self-bound structure at the following time U+i; 

(iii) the coalescence of the subhalo with the centre of the 
main halo (t en d). To compute this time, we adopt a dynam- 
ical friction time formul a just after the subhalo has been 
dissolved in episode (ii) |Binnev Sz Tremainelll987t ): 



tdS = Ofric(Sorb) 



Gm sub In A' 



(1) 



where afric(Oorb) encloses information on the subhalo orbit, 
Vvir is the virial velocity of the host, m su b and r su b are the 
mass and position of the subhalo relative to the host just be- 
fore dissolution, and In A = (1 -)- M v ir/m su b) is the Coulomb 
logarithm with M v i r th e virial mass of the host. W e take 
afric(e orb ) = 1.177? ' 78 (|Bovlan-Kolchin et all l200ll . where 
r\ = j/jc (E) is the orbital circularity of the subhalo relative 
to the halo centr^E The final cosmic time of the halo-halo 
central merger, i e nd, is the sum of tdiss+^df; we consider that 
t en d is a good approximation to the actual galaxy-galaxy 
merger epoch. 



2.2 Subhalo merger histories and central merger 
times 

Since our goal is to analyse the impact of mergers in the 
growth of the spheroidal component (bulges) of central 
galaxies, we extract the subhalo merger histories of the prin- 
cipal branches of a population of main subhaloes defined at 
z — 0. A main subhalo is the most massive structure within 
the hierarchy of subhaloes of a friend-of-friends (FOF) halo, 
and henceforth we refer to it as a distinct or main halo. We 
explicitly reject all subhaloes that coexist within a main halo 
at z — as being part of its merger history; although they 
will likely merge with the central object in the future, these 
subhaloes have no impact for the central bulge at z = 0. We 
also discard mergers within progenitors not associated with 
the main branch of the host halo. By following this proce- 
dure, our analysis is closer to the actual merger history of a 
population of central galaxies at z = than an alternative 
method based on FOF halo merger histories. 

We have randomly selected two samples of main 
haloes at z = 0, having 1347 and 1500 members with 
masses larger than 1.2 x 1O 12 M (10 3 particles) and 9.4 x 
10 10 M o (10 4 particles) for the MS-I and MS-II, respectively. 
Both samples (properly normalized to account for the frac- 
tional volume they cover relative to the whole simulation 
boxes) follow the mass function of the full halo population. 
For each halo in the samples, we obtain their subhalo merger 
history using the MS on-line databa sefl For details in th e 
construction of these merger trees see lSpringel et all (2005). 

A given merger event is denned by three epochs: 



1 http://gavo.mpa-garching.mpg.de/Millennium/ 



2.3 Galaxy occupation 

Although the merger histories of haloes may provide a basis 
for determining the morphology of their central galaxies, 
the impact of a merger does not depend directly on the 
halo mass ratio at the start of the merger but rather on the 
central dynamical masses (inner dark matter, gas and stars) 
that interact in the final stages of the galactic merger. 

To follow the stellar and gas mass assembly of galaxies 
inside the main haloes, as well as the processes that affect 
the gas and stars during me rgers, we use a semi-empirica l 
approach close to that one in Hopk ins et all (|2009bl . l2010ah . 
This approach yields stellar mass assembly histories that 
are consistent, by construction, with observational trends. 
In Appendix [A] we describe in detail this method and its 
implementation. In summary, for each present-day main halo 
in our MS samples we do the following: 

(i) extract the main branch of its merger tree; 

(ii) seed a central galaxy, the primary, at z se ed = 3 with 
stellar and gas masses given by the semi-empirical relations 
M»(M h ,z) and M 9 (M*,z) (Appendix |AT|) : 

(iii) identify :z s tart for each halo merger along the main 
branch (zstart ^ ^seed) and assign a galaxy to the in- 
falling halo, the secondary, according to M* (Mu, Zstart) and 
M 9 (M*,z start ); 

(iv) either assume that the stellar and gas content of the 
secondary does not evolve until final coalescence, or follow 
its evolution (the star formation, SF, and SF-driven outflow 
processes) by means of semi-analytic recipes (Appendix lASI) : 

2 The subhalo has specific angular momentum j and energy E, 
and jc(E) is the specific angular momentum of a circular orbit 
with the same energy and with a radius r c (E). 
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(v) compute the galaxy-galaxy (central halo) merging 
time, t cn d = idis + tdf, using Eq. (0 for tes; 

(vi) estimate the bulge and disc masses of the primary 
galaxy after coalescence at i en d, using physical recipes cali- 
brated by numerical simulations (Appendix I A2|l ; 

(vii) repeat the process for the central galaxy until reach- 
ing z — 0, taking care of each merger and updating at 
each z the properties of the central galaxy according to the 
M*(M h ,z) and M g (M,,z) relations. 

Note that in our scheme galaxies are initially pure discs, 
and that the only channel of bulge growth is through merg- 
ers. We begin testing our model by comparing the evolution 
of the galaxy merger rate with observational constraints. 



3 RESULTS 

3.1 Stellar merger fractions 

Fig.[T]shows /merg, the fraction of main halo progenitors that 
suffer at least one central stellar merger with /i„ > 0.1, rel- 
ative to the total number of main progenitors in a temporal 
bin of 1 Gyr, i.e., /merg is equivalent to the merger rate per 
halo and per Gyr. The blue, black, and red lines (circle, tri- 
angle and star symbols, respectively) are the median values 
for three different halo mass bins, denned at z = 0, centered 
in 3.3 x lO n M , 2.7 x 10 12 M© and 6.3 x 10 13 M Q , respec- 
tively; the error bars are the corresponding la scatters of 
the distributions. The mergers (including their stellar mass 
ratio) are defined at the time when the secondary galaxy fi- 
nally coalesces with the central one (z e nd)- In Fig.[T]we have 
actually used z S ced = 4, which is the maximum redshift used 
to infer the halo-to-stellar mass relation. 

The figure shows that the fraction of progenitors with 
on-going mergers grows with z, which is a trend inherited 
from the well known growth with z o f the halo-halo merge r 
rates (e.g., see right panel, Fig. 3 of iFakhouri et ail boiO). 
We can also see that /merg depends on Mh- galaxies in the 
most massive haloes have in average a higher fraction of 
mergers than those in the less massive ones; at the level of 
halo mergers, this dependence is quite weak. 

It is interesting to remark , as noted elsewhere 
|Stewart et al.1 120091 ; iHopkins et al J l2010al ) . that there is a 
dramatic change in the merger fraction as a function of mass 
and z depending in which mass ratio is used to define the 
merger. If we use fi, instead of /ihaio, there is an overall 
drop in the fraction of main progenitors with major merg- 
ers, and the dependence of /merg on mass increases, specially 
for the largest masses. This is mainly produced by the shape 
of the M*(Mh, z) relation: at masses below the knee of this 
relation, a /i,(z s tart) = 0.1 merger is actually related to a 
Mhaio(zstart) = 0.36 merger (because M h oc M° 44 ), which 
is more rare than a 1 : 10 halo merger; at masses above 
the knee, where Mh oc M, 2 , an opposite behavior is in prin- 
ciple expected. The stellar growth of the merging galaxies 
from Zstart until z e nd also contributes to the changes in the 
merger fraction. To give an impression of this, we show in 
Fig. [2]the median values of Phaio(^start)/M*( 2: cnd) as a func- 
tion of z for the mergers corresponding to Fig. [TJ in the case 
where the secondaries evolve according to a simple SAM (see 
Appendix lA3|) . The MS samples were divided into the same 
three present-day Mh bins as in Fig. [1] The error bars show 
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Figure 1. Relative fraction of main progenitors having at least a 
stellar major merger event with fi* > 0.1 as a function of redshift. 
The time interval between redshift bins is 1 Gyr. The fraction is 
relative to the total number of main progenitors in a given redshift 
bin. The colors (symbols) are for three different halo mass bins at 
z = according to the legend. The low-mass bin is for the MS-II 
sample, the other two are for the MS-I sample. The stellar major 
merger ratio and redshift are defined at the time of final galaxy 
coalescence. The secondary galaxies grow according to the SAM 
described in Appendix IA3I Stellar and gas masses for each halo 
are given by empirical relations described in App endix ) A 1J The 
shaded region is a compilation of observations bv lHopkins et all 
teoiod lbl) for £t„ > 0.1. The magenta arrows mark the change on 
/merg for the intermediate mass bin if the major merger threshold 
is increased to 0.25. 



the la scatters. On average, a /ihaio(zstart)=0.1 halo merger 
corresponds to a ^i*(z C nd) ~ 0.01 stellar merger for low- 
mass galaxies. Interestingly, notice that in our model of sec- 
ondary evolution, the /i* ratios are roughly the same whether 
they are measured at the beginning (z sta rt, dashed lines) or 
at end of the merger (zend, s °hd lines). This means that 
in average, the secondary increases its stellar mass (semi- 
analytically calculated) roughly by the same amount as the 
primary (semi-empirically assigned) during the period from 
infall to coalescence. 

In the case of the baryonic (stars + cold gas) merger- 
mass ratio, the fraction of haloes having mergers with /ibar > 
0.1 increases relative to that based on stellar major merg- 
ers. This is because the fraction of cold gas is significantly 
high for galaxies with low stellar masses, specially at high 
redshifts. A significant fraction of mergers with /j,,< 0.1 
are actually major baryonic mergers just because the sec- 
ondaries have larger cold gas fractions than the primaries. 
Morphology-based techniques for estimating merger rates 
are more sensitive to baryonic than to stellar mergers, and 
indeed they find higher mer ger rates than with pair-based 
techniques l|Lotz et al.ll201ll ). 

The shaded area in Fig. [1] encomp asses the ob- 
served major-merger fractions compiled by Hopk ins et al.l 
(|2010al lbh. Most of the observational results do not have a 
well-defined mass selection and are based on different merger 
identification criteria. Very roughly, the stellar mass of the 
primary galaxies covers a range 10 10 — 2 x 10 11 Mq at low 
redshift, with a larger minimum mass at high z. For mea- 
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Figure 2. Median values and la scatter of the distribution of 
halo- to-stellar merger mass ratios, /^halo (zstart ) /m* (z C nd ) > f° r the 
mergers shown in Fig. [l] with the same three halo mass bins 
(solid curves from top to bottom are for lower to higher masses, 
respectively) and for the case of evolution in the secondaries. The 
dashed lines are the same but with /i* defined at the beginning 
of the halo merger. 



surements based on pair samples (pre-merger), ^t* > 0.25, 
while for measurements based on morphology samples (post- 
merger), /i» > 0.1. The lower (upper) bound of the shaded 
region is dominated by the pair (morphology) samples, and 
therefore reflects merger fractions with /i* > 0.25 (/z* > 0.1). 

Fig. [I] shows that there is a reasonable agreement be- 
tween predictions and observation^]. If we limit the mergers 
to (U*(z e nd)> 0.25, then the fraction falls below the lower 
limit of observational estimates for z < 1 (magenta arrows 
for the intermediate mass bin where most galaxies have 
M*(z = 0)~ 1 - 5 x 1O 1O M ). In a recent study based 
on the 30-band photometric catalogue in COSMOS, com- 
plemented with the spectroscopy of the z COSMOS survey 
to define close pairs. [L6pez-Saniuan et ail l|2012f i report the 
merger fraction evolution (0.2 < z < 1) of massive galax- 
ies, M* > W^Mq (corresponding to M h > 7 x 1O 12 M ), 
for fi, > 0.1. The predicted merger fraction evolution is in 
reasonable agreement with these observations. 



3.2 Merger-induced bulge growth and its 
dependence on satellite evolution 

After showing that the ACDM-based galaxy merger frac- 
tions (rates) are roughly consistent with observations, we 
concentrate now in the growth of bulges during mergers 
through the following processes (see Appendix I A2|) : 

(a) the total acquisition of stars from the secondary, 

3 We note that while our predictions refer to evolutionary tracks 
of individual galaxies defined at z = 0, observations refer to sam- 
ples above a constant M„ at all epochs. For very massive galaxies, 
M, changes little since z ~ 2, but for intermediate and low-mass 
galaxies, M* is smaller at higher redshifts. Hence, the prediction 
for the merger fraction above a certain M* threshold, constant at 
all epochs, would be curves slightly steeper than the ones shown 
in Fig[l] particularly for the intermediate and low-mass bins. 



(b) the violent relaxation and transport to the centre of a 
fraction of stars in the primary disc, and 

(c) the newly formed stars in central starbursts produced 
by a fraction of the gas from the merging galaxies (Eq. A2). 

To calculate the amount of stars added to the bulge 
through processes (a)-(c), we estimate the dynamical (stars, 
gas, inner dark halo) masses of the merging galaxies. As 
mentioned above, the empirical M*(Mh,z) and M 9 (M», z) 
relations are used for this whenever is possible. These re- 
lations hold for the primary at any point before final coa- 
lescence and shortly afterwards. For the secondary however, 
they are only valid before the start of the halo merger. As 
the secondary spirals inwards towards the centre, different 
physical processes alter its stellar and gaseous contents. To 
study how important some of these processes are, we first 
analyse the case where the gas and stellar masses in the sec- 
ondary are "frozen" at the values they had at the time the 
halo merger started (§ £|3.2.1[1 . and we later move to the more 
realistic case of stellar and gaseous evolution in the satellite 

(§33321. 

In order to use both MS data-sets without introduc- 
ing a bias due to numerical resolution, we consider only 
those mergers with fthaLo ^ 0.1 in what follows. At this level, 
both simulations capture the majority of the halo mass con- 
tributed by mergers and statistically match each other in 
their overlapping mass range. This could raise the concern 
that by missing all mergers with /^haio < 0.1, the growth 
of bulges might be affected significantly. However, due to 
the shape of the M*(Mfe, z) relation, halo-halo mergers with 
A*haio = 0.1, imply galaxy-galaxy mergers with /i* <C 0.1 
at low masses and /i» ~ 0.1 at large masses (see Fig. [2}. 
Thus, it is only for the most massive systems that some rel- 
atively large stellar mergers are being excluded, which are 
in any case not very relevant since the bulk of the bulge 
is formed during the la rgest merger events (see e.g. Fig. 8 
of iHopkins et al]|2010a[ ). For massive haloes, we have nev- 
ertheless checked the role of including those events with 
Mhaio < 0.1, only possible in the MS-II, and found that the 
results change very little. 



3.2.1 Case 1: non-evolving secondary 

The upper left panel of Fig. [3] shows B/T as a function of 
M*(z = 0) for the central galaxies seeded into the MS-I 
and MS-II haloes (circle and star symbols, dotted and full- 
line shaded regions, respectively) in the case where secon- 
daries have, at coalescence, the same M* and M g they had 
when their haloes became subhaloes (extreme SF quench- 
ing). Since mergers are more frequent in massive galaxies 
and since the gas fraction is a decreasing function of M*, 
it is expected that massive galaxies will be more prone to 
develop large bulges. Several studies in the past have shown 
the beh avior depicted in Fig. [5] for a diverse class of mod - 
els (e.g. IHopkins et al.ll2009bl . l2010al ; I5e Lucia et ai1l201ll ). 
Below a few times 10 10 Mq, most galaxies have B/T ^ 0.2. 
This is partially an effect of the high gas fractions of low- 
mass galaxies, which make the bulge growth less efficient. 

The middle panels of Fig.[3]show the contribution of the 
three merger-driven processes of bulge growth mentioned 
above, as a function of M*. The bulge mass fraction as- 
sembled through processes (a), (b), and (c) is shown with 
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Figure 3. Le/i panels: Bulge-to-total stellar mass ratio as a function of M„(z = 0). Middle panels: Fraction of the final stellar bulge 
mass assembled by the three mechanisms of bulge growth in a merger: stars formed in starbursts (blue), stars transferred from the 
primary disc (red), and stars from the secondary (orange). Right panels: Comparison between the fraction of stars in the primary, 
that are transferred to the bulge and the fraction of stars that are added via the merger of the secondary, ff rl ; the one to one relation 
between these fractions is shown as a dotted line. These fractions are relative to the mass of the primary just before final coalescence. 
In all panels, the symbols and solid lines are the median values and the shaded areas contain the ±<r regions of the distribution. Circles 
(stars) and dotted (full) lines patterns are for the MS-I (MS-II). In the right panels, the sample of galaxies was divided in two mass 
ranges: M* < 1O 1O M0 (MS-II data only) and M* > 10 1:l Mq (MS-I data only). Galaxies were seeded at z scc( j = 3 and only mergers with 
Mhalo > 0.1 have been included. The upper row is for the case where the fraction of gas and stars in the secondary is the same at the 
time of final coalescence as it was when the halo merger started, whereas the lower row is for a model where the gas and stars in the 
secondary evolve through SF and SN feedback after the halo merger starts (Appendix IA3II . Observational data for two galaxy samples 
with M* ^ 1O 1O M0 from G09 and FD11 has been added to the left panel with red circles and blue squares, respectively. The median of 
the distribution for the galaxy samples are shown with solid lines. 



orange, red and blue, respectively. In the right panels of 
Fig. [3l we compare the fraction of stars violently relaxed in 
the primary (which is a function of the effective dynamical 
mass ratio, fi e ff, see Eq, IAl|) to the fraction of stars acquired 
from the secondaries (which is equal to the stellar mass ratio 
Ht) for each merger. We have divided the events according 
to the mass of their descendant at z = 0: massive galaxies 
with M* > 1O 11 M0 (circles and dotted region) and low-mass 
galaxies with M* < 10 10 M Q (stars and dashed region). 

In the case of no satellite evolution, the bulges of 
M, < 10 11 M e galaxies grow overwhelmingly through stars 
transferred from the primary rather than by stars acquired 
from the secondary with fi, < 0.1 (i.e., the stellar mergers 
are essentially always minor). This is mainly because the 
secondaries are gas dominated in most of the cases. For the 
same reason, we see in the middle panel (upper row) that the 
contribution from starbursts (blue) for galaxies of this mass 
is also more important than the one from stars in the sec- 
ondaries (orange). Nevertheless, starbursts never contribute 



more than ~ 10% to the present-day mass bulge, and they 
are essentially negligible for massive galaxies due to their 
very low gas content. 

For galaxies with M* > 2 x 10 11 M Q , it is possible to 
have mergers with fi t > 0.1, and when these mergers reach 

0.2, then the dominant channel of bulge growth is stars 
directly added from the secondary. Thus, the difference be- 
tween galaxies of different mass in the right panel of Fig. [3] 
explains clearly the trends seen in the middle panel be- 
tween the red and orange distributions. This is ultimately 
connected to the shape of the M*(Mh) relation. 

3.2.2 Case 2: evolving secondary 

The effect of including SF and SN feedback in the evolution 
of the secondary in its travel to the centre of the primary 
(Appendix IA3|) is explored in the lower panels of Fig. [3] 
These processes augment the stellar mass of the secondary 
and reduce its gas fraction. They are of course more rele- 



© 0000 RAS, MNRAS 000, 000-000 



The growth of galactic bulges through mergers in ACDM haloes revisited 7 



vant when the supply of gas at the start of the merger is 
high and thus, they considerably affect low-mass galaxies 
and are nearly irrelevant for massive galaxies. In what fol- 
lows we take this model as our fiducial one. Note that this 
case also implies some level of SF quenching since once a 
galaxy becomes satellite, there is no newly accreted gas. 

The B/T ratio of low-mass galaxies is on average 
slightly higher in this case. Continuous SF increases the stel- 
lar mass of the secondary at the time of coalescence relative 
to the start of the merger, making the contribution from 
stars in the secondaries to the bulge growth much more sig- 
nificant (orange region, middle panels of Fig. [3]). Galaxy 
outflows from SN feedback attenuate slightly this effect by 
removing gas from low-mass secondaries. On the other hand, 
the other two mechanisms (b) and (c), are diminished be- 
cause /ieff ~ Mbar is reduced due to the loss of gas (compare 
right panels of Fig. [3]). 



Bulge stars formed in merger-induced starbursts con- 
tribute on average less than ~ 5% to the total masses of 
present-day bulges, the largest contribution being in galax- 
ies with masses Af* « 2 x I0 10 Af©. H09a and lHopkins et all 
(2010a) report a larger contribution, partic ularly for low- 
mass g alaxies, reaching ~ 40% (see Fig. 14 of lHopkins et al.l 
2010a ). This d ifferenc e is partly caused by the model in 
Hopkins et al.l l|2010al ) being closer to our previous model 
where the gas fractions are not allowed to decrease while the 
satellite merges. Moreover, we also speculate that an addi- 
tional cause of the discrepancy is the different treatment of 
the extrapolation of the gas fraction f g = M g /M* to lower 
masses and higher redshifts: we have put an upper limit on 
f g based on the maximum observed value in th e compila- 
tion of o bservations give n bvlStewart et all (|2009l ): f B < 100, 
whereas iHopkins et al.l |2010al ) uses a less restrictive limit 
(see Appendix IA1|) . We found that the latter case indeed 
creates a larger contribution from starburs ts in low-mass 
;alaxie s, although not to the level reported in lHopkins et al.l 
2010al ). It also creates more destructive mergers (more mas- 
sive bulges), simply because there is more SF prior to the 
final coalescence. Because of this, it is hard to increase the 
contribution from starbursts while at the same time keeping 
the B/T distribution consistent with observations. 



3.2.3 Comparison with observations: B/T demographics 

We first remark that the bulge/disc decomposition 
of observed galaxies and the characterization of a 
spheroid as clas sical or p seudo - bulge are very d i fficult 



In spite that the contribution of stars from the secon- 
daries is increased compared to the case of no satellite evo- 
lution, the contribution from stars violently relaxed in the 
primary continues to be dominant for Af* ;S 10 11 Mq. For 
Af* » 10 9 , 3 x 10 10 , and 10 11 Af®, on average60%, 55%, and 
45% of the final bulge mass was accreted through this chan- 
nel, respectively. This implies that the smaller the galaxy, 
the more their bulges share properties of their stellar popu- 
lation with their discs. Notice that this is an important result 
since it suggests that bulge growth due to secular evolution 
and that induced by mergers cannot be distinguished only 
by looking at the similarity of the stellar populations of the 
disc and the bulge, as is sometimes assumed. There are how- 
ever some small galaxies with bulges mostly assembled by 
stars from the secondaries; we can see that this happens 
when the mergers with /i* > 0.25 dominate (right panel of 
Fig. O lower row) . 



Allen et al.l 2006|; Laurikaincn et al. 2007 


; Fisher & Drorvl 


20081; iTasca & White! l201ll; ISimard et al.l 


20111). There are 



only a few studies of bulge/disc decomposition applied to 
volume limited samples that can be used to obtain fair statis- 
tics on the B/T distribution as a function of Af*. We com- 
pare our predictions with a couple of these studies in the 
following. 

In the lower left panel of Fig. [3] we have included two 
observational samples. The first (shown in red) is a volume- 
limited sample of ~ 1000 local galaxies from the Sloan Dig- 
ital Sky Survey (SDSS) with Af* > 1O 1O M taken from 
iGadottil (|2009l ) (hereafter G09). Since we are modeling cen- 
tral galaxies only, we have removed satellites from this sam- 
ple by cross-correlating the original sample with the S DSS 
galaxy group catalog constructed by I Yang et al.l (|2007l ). We 
find that 25% of the original sample are satellites, 64% cen- 
trals, and the rest are unidentified. We keep the latter two for 
comparison with our mock sample (we note however that we 
did not find a significant difference between the B/T distri- 
bution of central/unidentified galaxies and that of satellites: 
the latter is just slightly more peaked at B/T ~ 0.2 — 0.3 
than the former). The median values of B/T for the sample 
of central galaxies (solid red circles) increase with mass more 
steeply than our predictions, with median values larger by 
factors up to 1.5, although the scatter in the observational 
sample is very large. In particular, there is a significant frac- 
tion of observed galaxies with B/T = 1. 

The second sample (shown in blue) is from FD11, which 
is a volume- limited sample of galaxies within the local 11 
Mpc volume. We have selected only those galaxies with 
Af, 10 9 M© (a total of 99 galaxies). Unfortunately, for this 
sample we do not know which galaxies are central/satellites. 
At masses 9.8 <log(Af*/Af©)< 10.8, the observed median 
values of B/T are close to those of our predictions, well 
within the statistical scatter. However, at lower masses, the 
median B/T ratios of our mock galaxies tend to be larger 
than those of the FD11 sample. 

In the top panels of Fig.|4]we show the B/T distribution 
for galaxies with Af* above a certain threshold (1O 9 M0 and 
10 10 Mq for the left and right panels, respectively) for a com- 
bination of the MS-I and MS-II samples. To obtain this com- 
bined MS sample, we used the latter up to Mh = lO 12 Af0 
and the former for haloes with higher masses but with a 
number count renormalized to match the mass function 
of the MS-II sample in the mass bin 1 - 3.6 x 10 12 M Q , 
where both samples have neither a completeness nor a low 
number statistics problem. We have checked that the re- 
sults do not change significantly if the samples are taken 
separately in their respective mass range of completeness. 
The statistical error bars shown in th e figure are Poisso - 
nian, employing the definition given bv lLukic etal] <|2007l ): 
a± — y/N~+T/4 ± 1/2, where N is the number of main 
haloes in a given bin. We also show the observational sam- 
ples from FD11 (blue), G09 (red) and W+09 (magenta). 
These panels show what was apparent in the left panels 
of Fig. [3] the B/T distribution strongly depends on the 
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Figure 4. Top left: B/T distribution for galaxies with M t ^ 10 9 Mq in the combined MS sample (black histograms) and in the 
observational sample from FD11 (blue symbols). Top right: The same as top left but for galaxies with M* > 10 10 Mq; the galaxy samples 
from G09 and W+09 have been added to the figure with red (star) and magenta (triangle) symbols, respectively. The latter is compared 
with the models only for B/T < 0.5. Errors in the number counts are Poissonian and are marked with bars for the observational data 
and with shaded regions for the model predictions. Bottom left: The combined MS sample with these different mass thresholds is divided 
in galaxies with a "classical bulge" , defined as those where the stars transferred from the primary contribute less than 50% to the final 
bulge mass; the remaining galaxies are called "pseudo-bulges" by extension. The panel shows the B/T distribution for classical bulges 
with a solid (dashed) line for M t > 10 9 (10 10 )Mq. The observational samples of G09 (red stars), FD11 (blue circles) and W+09 (magenta 
triangles) are also shown. Bottom right: Fraction of "bulgeless" galaxies as a function of M*, with "bulgeless" defined as those galaxies 
with B/T < 0.1 (red), and these plus all pseudo-bulges with B/T > 0.1 (black) 



The simulation data is shown with histograms and the observational data with stars and circles for G09 and FD11, respectively. 



mass cut-off of the sample. For representative samples with a 
M* threshold of 1O 9 M , which are dominated by low-mass 
galaxies, the fraction of galaxies with B/T < 0.1 is 20% 
higher than for samples with a threshold of 10 10 Mq. 

The B/T distribution of our mock galaxies roughly 
agrees with the one from the FD11 sample; although the 
latter has a higher fraction of galaxies with B/T < 0.1 
and slightly lower with 0.6 < B/T < 0.9. We recall, how- 
ever, that the FD11 sample is very local and hence, mas- 
sive galaxies (with typically high B/T values) are under- 
sampled. These massive galaxies are represented better in 
the G09 sample, which shows a flatter B/T distribution 
than the one from FD11, and compared to our results, it 
has a smaller fraction of galaxies with B/T < 0.2 and a 
larger fraction of galaxies with B/T > 0.6. The sample of 



W+09 has a morphological selection criterion that excludes 
elliptical galaxies, and this is why we have restricted the 
comparison to B/T ^ 0.5; for lower values of B/T, this 
sample has a similar distribution to the one from FD11, and 
thus, agrees with our predictions. If we assume that bulge- 
dominated galaxies (B/T ^ 0.5) correspond to E/SO mor- 
phological types, then the fraction of such types predicted 
by our model is « 16% for M» > 1O 1O M . This is slightly 
lower than what was found observationally in lBaillard et al.l 
Hon]), where » 20 - 25% of galaxies are of E/SO type at 
z ~ 0.015 (for a sample with a completeness above 80% for 
M, > 1.5 x 10 10 M Q ). Our fraction could be larger if one 
considers that some SO gal axies actually have B/T < 0.5 
l|Kormendv fc Bendeulioia ). 

We note that our model is not able to produce galax- 
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ies with B/T > 0.9. For the most massive galaxies, this 
could be because the primordial galaxies seeded at z SC ed = 3 
are assumed to be discs, which implies that some fraction 
of these initial discs survives even after major mergers. In 
fact, several pieces of evidence suggest that the most mas- 
sive galaxies (ellipticals with M* > 5 x 10 11 A/0) formed 
very early (z > 2) and have not grown in mass by in situ SF 
since then (e.g.. iThomas et al . I [20051 ; I Collins et alj|200gh . If 
we were to assume that the most massive of our modelled 
galaxies where spheroids instead of discs at z acc d, we would 
have a few B/T ~ 1 massive galaxies. 

We also remark that in our scheme a disc may re- 
main even after a merger with a very large mass ratio, con- 
trary to some S AMs that assume that in these cases all the 
baryons end in a st ellar spheroid (e.g.. IParrv et al.ll2009l ; 
|Pe Lucia et alJboill ). It should also be noted that (i) the 
semi-empirical M,(Mh) and Af ff (Af*) relations refer to av- 
erage trends where additional dependencies on environment 
are not considered, and (ii) we do not follow the morpho- 
logical ev olution of satellite gala xies. According to the SAM 
results in Lucia et al.l (|201 lh . the latter is not relevant 
to obtain B/T > 0.9 galaxies. Instead, the possibility of 
disc conversi on into a spheroid due to intrinsic disc instabil- 
ities (see e.g. IParrv et alj|2009h . which is something we did 
not take into account in our modelling, contributes to the 
difference with our predictions. While intrinsic secular evo- 
lution is a reasonable mechanism to produce pseudo-bulges 
in dynamically cold discs, it is a matter of debate if it could 
transform most of a massive stellar disc into an elliptical 
galaxy or destroy a small disc embedded into a dynamically 
hot spheroid. 

On the observational side, it is important to notice 
that in many cases of bulge/disc decomposition in ellipti- 
cal galaxies, adding a small disc does not improve the sta- 
tistical significance of a pure bulge model, and a value of 
B/T — 1 is simply assigned. However, several studies pro- 
vide evid ence to argue that most ellipticals cont ain small 
discs (e.g. iFerrarese et al.lfl99ll ; lLauer et al.ll 20051 '). 

Fig. [5] is another way of showing the predictions of our 
model. Here we plot the fractions of galaxies in different 
M* bins divided according to their B/T ratios: B/T < 0.1 
(filled circles, red line), 0.1 ^ B/T < 0.5 (stars, blue line), 
and B/T ^ 0.5 (squares, black line). Roughly, these three 
groups can be associated to irregulars (types later than Sc), 
intermediate disc-dominated (Sa to Sc types), and bulge- 
dominated (E/S0) galaxies, respectively. At M* < 10 10 M Q , 
w 65% of our model galaxies have B/T < 0.1, at M* f» 
3 - 8 x 10 10 M Q more than 50% have 0.1 B/T < 0.5, 
while at M» > 10 11 M©, galaxies with B /T ^ 0.5 dominate. 
These predictions are roughly consistent with estimates of 
the local morphological mix. The predicted mass fractions of 
stars contained in discs (Af* > 10 9 Mq) is 57% and the rest 
is in spheroids, consistent with observations. The fraction of 
stars in galaxies with B/T ^ 0.1 is « 15%, while the fraction 
of stars (both in the disc and the spheroid) in galaxies with 
B/T > 0.5 (spheroid-dominated) is « 58%. 

3.2.4 Comparison with observations: bulge composition 

We study the fraction of galaxies whose bulges were assem- 
bled mostly by stars from the secondaries (contributing to 
the bulge mass fraction by more than 50%). This fraction as 



a function of B/T is shown in the bottom left panel of Fig. 2] 
(solid and dashed lines for A/* > 1O 9 M and > 1O 1O M , 
respectively). Because this channel of bulge growth is dom- 
inant when stellar major mergers dominate the bulge mass 
assembly (see right panels of Fig.[3|, we will nominally refer 
to these bulges as classical-like (CL). On the contrary, those 
bulges where more than 50% of their stars come from the pri- 
mary will be nominally defined as pseudo-like (PL) bulges. 
Notice that this division is well motivated: classical (pseudo) 
bulges are thought to be formed by major mergers (disc sec- 
ular evolution/ minor mergers). Thus, one expects that clas- 
sical (pseudo) bulges are dominated by stars acquired from 
the secondaries (primary discs). Moreover, the kinematics of 
the bulge stars, often regarded as one of the key discrimi- 
nants between pseudo and classical bulges, are expected to 
be consistent with this division: (i) bulges whose stars come 
predominantly from the secondaries will have more random 
orbital orientations (not necessarily aligned with the disc) 
and a bias towards radial anisotropy, (ii) bulges where most 
of the stars were scattered by instabilities resonantly excited 
by mergers (or formed "in situ" in the starburst) will have 
the bulge preferentially aligned with the disc, with more ro- 
tation, and tangentially biased orbits. 

We find that most of the bulges in our mock galaxies 
with B/T > 0.25 are CL, while the ones of small galaxies are 
typically PL. For B/T < 0.1, only < 25% of the bulges are 
CL. The red stars, blue circles and magenta triangles on this 
figure correspond to the observational samples of G09, FD11 
and W+09, respectively, where the bulge class is determined 
from the light profile (photometry). Our prediction roughly 
agrees with these observational samples, although it slightly 
under(over)-predicts the abundance of classical bulges in the 
intermediate (low) B/T range. 

In the bottom right panel of Fig. [4] we show the fraction 
of galaxies with B/T < 0.1 (which we call "bulgeless", red 
histogram) as a function of M* . For the black histogram we 
add to this fraction those galaxies with B/T ^0.1 that have 
PL bulges. The circles and stars correspond to the respec- 
tive fractions for the samples of FD11 and G09, respectively. 
The deficit of CL bulges for intermediate B/T values shown 
in the lower left panel appears in this figure as an overpre- 
diction of PL bulges for galaxies with M* ~ 5 x 10 10 Mq by 
a factor of ~ 1.5. 

Our model seems to underpredict the fraction of bulge- 
less galaxies by less than 15% at M» < 10 10 Mq. At larger 
masses, the agreement is good although there is an excess of 
bulgeless galaxies, which actually grows to more than a fac- 
tor of ~ 2, once we compare with the G09 sample; compared 
to the W+09 sample however (M» > 10 10 Mq disc galax- 
ies), our predicted fraction is slightly lower. Thus, it seems 
that the potential problem of too few bulgeless galaxies is 
related only to galaxies with Af* < 10 10 Mq . We recall that 
the FD11 sample contains both central and satellite galaxies, 
while our results refer only to the former. Because centrals 
have a higher probability to suffer mergers than satellites, 
with the difference being larger for smaller A/*, we expect 
that bulgeless galaxies are more frequen t in low-mas s satel - 
lites than in centrals. Indeed, using the lWang et al.l (|2008T ) 
SAM and the H09a procedure to estimate bulge growth, 
iFontanot et al.l (|201lT ) found that, at a given mass, low-mass 
satellites are more likely to be bulgeless than centrals. Un- 
fortunately, we do not know the fraction of satellites in the 
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Figure 5. Predictions of the relative fraction of galaxies as a 
function of M* according to the final bulge-to-total mass ratio: 
B/T < 0.1 (red circles), 0.1 < B/T < 0.5 (blue stars) and B/T > 
0.5 (black squares). 

FD11 sample, but if we could exclude it, it is likely that the 
bulgeless fraction would be lower than the one shown in the 
bottom right panel of Fig. [3] 

The predicted fraction of bulgeless galaxies could also 
decrease if one takes i nto account the effect of intrinsic disc 
instabilities, although lFontanot et al.l (|201ll ) found that this 
is only a minor effect for low-mass central galaxies. Despite 
that this effect is more relevant for massive galaxies, most 
of these already have large B/T (> 0.1) ratios due to major 
mergers. 



4 SUMMARY AND CONCLUSIONS 

We present a scenario of bulge growth in central galax- 
ies based on: (%) the mass aggregation and merger histo- 
ries of ACDM distinct haloes randomly selected from both 
Millennium simulations (§2), (it) the use of ^-dependent 
empirical relations, M,(Mh) and M g (M*), to seed galax- 
ies into the growing (sub)haloes (§2.3, Appendix Al), and 
(Hi) the implementation of a physically motivated model 
to calculate the bulge stellar mass after a galaxy-galaxy 
merger (Appendix A2). By construction, the stellar mass 
assembly of our mock galaxies follows on average the 
general downsizing trends in ferred from observations (see 
iFirmani fc Avila-Reesell2010r ). 

In our approach, bulges grow through galaxy mergers 
characterized by the dynamical mass ratios of the merging 
galaxies at the coalescence time. The latter is estimated fol- 
lowing the orbit of the subhalo until its number of parti- 
cles falls below the resolution of the simulations, and later 
applying a tested approximation of the dynamical friction 
time (Eq. [1]). The stellar and gas masses of the secondaries 
at the moment of coalescence are either assumed to be the 
same they had when their haloes became subhaloes (extreme 
satellite quenching) or their further evolution until coales- 
cence is estimated using a simple model of SF and SN-driven 
outflows (Appendix IA3|) . under the assumption that there 
is no further gas accretion. The bulge stellar mass increases 



after a merger by three contributions: (1) all the stars from 
the secondary, (2) a fraction of stars transferred from the 
primary disc, and (3) stars formed in a central starburst 
triggered by the condensation of cold gas from both merg- 
ing galaxies. While the first contribution is expected to be 
related to a classical bulge, the other two are associated to 
a pseudobulge, i.e. a bulge that shares properties with the 
disc. 

Our goal has been to revisit and explore in detail the 
predictive power of the ACDM paradigm to explain, through 
the halo-halo merger histories mapped in a non-trivial way 
to galaxy-galaxy merger histories, the properties and demo- 
graphics of the bulges of central galaxies with M* ^ 10 9 
Mq. We highlight the following results and conclusions: 

• Mock galaxies with M» < 2 x 10 10 Mq have on average 
B/T < 0.2, while for more massive galaxies, B/T increases 
rapidly with M*. This is mainly because the stellar major 
merger rate history increases with mass. For M* < 10 10 
Mq, nearly 60% of galaxies have B/T < 0.1, while for 
M» > 10 11 Mq galaxies with B/T > 0.5 dominate. The 
driver of bulge growth in massive/high-_B/T galaxies are ma- 
jor/intermediate (/i* > 0.1) stellar mergers. Instead, for low- 
mass/low-B/T galaxies, bulges grow mainly through mi- 
nor/minuscule stellar mergers (/i* < 0.1), although the cor- 
responding dynamical merger ratios are larger, /i c ff > /i*. 
The latter implies that the bulges of < 10 11 Mq galaxies 
have grown mostly from stars of the primary disc. 

• Bulges are composite, growing in several episodes 
through concomitant channels of stellar mass acquisition: 
from the secondaries, from the primary disc, and from star- 
bursts. The first (second) of these contributions increases 
(decreases) monotonically with M, . For M» « 10 9 , 3 x 10 10 , 
and 10 11 M Q , approximately 60%, 55%, and 45% of the 
bulge mass was assembled by stars from the primary, re- 
spectively. For M* > 10 Mq, the stars acquired from the 
secondaries contribute > 50% to the bulge mass. Bulge stars 
formed by central starbursts contribute < 5% to the bulge 
mass on average. Although intrinsic disc instabilities were 
not taken into account as a mechanism of (pseudo)bulge 
formation, the merger-induced disc instabilities produce a 
similar effect regarding bulge mass acquisition. Since this 
mechanism dominates at low M*, the stars of the bulges of 
small galaxies are predicted to have similar properties, in- 
cluding their kinematics, than the stars of their discs. We 
can nominally consider that when stars from the secondaries 
contribute more (less) than 50% to the final bulge mass, 
then the bulge is classical-like (pseudo-like). The fraction of 
classical-like bulges is > 50% for galaxies with B/T > 0.25. 
Galaxies with smaller B/T values have mostly pseudo-like 
bulges. 

• The evolution of the satellites from the time of infall 
until final coalescence influences the B/T ratio of the central 
galaxy and significantly affects the ratio of secondaries-to- 
primary stars in the remnant bulge. In the extreme case 
of total SF quenching in the satellites, low-mass galaxies 
end up with smaller B/T ratios than in the case where 
SF/feedback is allowed. In the former case, the contribu- 
tion of stars from the secondaries to the growth of the bulge 
falls to < 10%, resulting in a negligible fraction of classical- 
like bulges for galaxies with M, < 10 11 Mq, which is in 
conflict with observations. On the contrary, if satellites keep 
and transform their gas into stars too efficiently, then the 
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mergers would be very destructive and the central galaxies 
would end up with B/T ratios that are too high. An in- 
termediate case between the previous two extreme cases is 
the standard SAM we have used to follow SF/feedback in 
the satellites (see Appendix IA3|) ; it assumes that their hot 
gas is completely stripped as they merge with the hoslQ. For 
this case, it turns out that on average, the stellar mass of 
a satellite of a given mass grows in a similar way as that 
of a central galaxy of the same mass, with the latter being 
determined by the evolution of the M*(Mh) relation. Since 
this case is quite successful, as we summarize below, the lat- 
ter finding is also very useful as a reference to infer merger 
rates from observations based on pair samples at relatively 
large separations. 

Are our results in agreement with observations of 
bulge/disc demographics? Firstly, we checked that the pre- 
dicted stellar merger rates as a function of mass and z are 
actually consistent with observational measurements (Fig. 
[TJ. Then, we compared our results with the observational 
samples of FD11, G09, and W+09 regarding the B/T de- 
pendence on M* (Fig. [3]), the B/T distribution and the frac- 
tions of galaxies with a given B/T and M* (including the 
fractions of classical- and pseudo-like bulges; Fig. De- 
spite the large observational uncertainties, the overall con- 
sistency between our ACDM-based predictions and observa- 
tions is remarkable. In particular, it is quite relevant that 
a large fraction of the pseudo-like bulge population can be 
explained in a merger-induced scenario without introducing 
pure intrinsic secular instabilities on the disc. We also find 
an agreement in the fractions of galaxies within a given B/T 
interval, which can be roughly associated to different mor- 
phological types (the morphological mix, Fig. O, as well as 
in the fractions of stars in the disc and bulge components. 
While overall our prediction for the B/T distribution is con- 
sistent with observations, in detail we have detected some 
potential disagreements: 

-There are no mock galaxies with B/T ~ 1, while ob- 
servers assign B /T — 1 to a small fraction of (discless) galax- 
ies. We argue that this apparent disagreement rather than 
a failure of the model could be a prediction: even the giant 
ellipticals should have small hidden discs. 

-At the mass range 10 9 < M,/M Q < 10 10 , there arc 
15% fewer mock galaxies with B/T < 0.1 than observed 
in the local 11 Mpc volume sample of FD11. The latter, how- 
ever, is likely to contain a non-negligible number of satellites 
that increase the fraction of bulgeless galaxies (satellites are 
expected to have smaller B/T ratios than centrals of the 
same mass). On the other hand, the fraction of low-mass 
B/T < 0.1 mock galaxies can easily be increased if the low- 
mass end of the M, [Mh , z) relation is such that at higher 
redshifts the values of M* for a given Mh are smaller than 
those used here (see Appendix B). 

We conclude that the implementation of a reasonable 
semi-empirical model of galaxy occupation into growing 
ACDM haloes is able to predict the present-day demograph- 
ics and mass dependence of the galaxy B/T ratios. At the 



4 If environmental processes such as ram pressure stripping not 
only remove all the hot gas from the satellites but also their cold 
gas discs, then star formation will be further suppressed and it 
would be more difficult to reproduce the observations. 



basis of our results are the ACDM halo merger rates, the 
merger-driven bulge growth model, and the way galaxies as- 
semble their gas and stars, with a downsizing trend and a 
bell shaped (low) stellar formation efficiency. 

Although these results might be interpreted as a pos- 
itive test of the ACDM model, we remark that we have 
only included merger-driven mechanisms of bulge growth. 
Discs instabilities such as bars and spiral arms also lead to 
secular processes of pseudobulge growth, in particular, the 
" peanut-shaped" /boxy b ulges are thought to be a byprod- 
uct of bar formation (e.g. iKuiiken fc Me rrificld 199^). There 
are other two possible mechanisms of bulge formation: (i) 
fragmentation of the gas-rich disc into clumps that migrate 
towards the cente r forming a large bulge at high redshifts 
jDekel et al.ll2009l . a lthough stellar feedb ack might suppress 
this mechanism, see lHopkins et al]|201 j ). and (ii) misalign- 
ment of the angular momentum of th e disc with that of the 
newly-accreted gas (Sale s et al.ll20lih . All these non- merger 
processes would probably increase the B/T ratios of our 
mock galaxies, decreasing the population of bulgeless galax- 
ies. Since our predictions are already in marginal agreement 
with observations of low-mass galaxies, the addition of these 
mechanisms might pose a challenge for the ACDM model. 
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APPENDIX A: SEMI-EMPIRICAL MODEL OF 
GALAXY OCCUPATION 

Al Stellar and gas mass relations 

Along the evolution of each halo we assign stellar masses 
by using semi-empirically determined average stellar-to- 
halo mass relations at different epochs. By means of the 
abundance matching technique, iBehroozi et al.l (|2010t ) de- 
termined the M,-Mh relations for central galaxies in the 

< 2 < 1 range and then extended these relations to 

1 < z < 4. These authors presented two different fitting for- 
mulae to their results, one for each redshift range. Because 
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these formulae ar e disjoint at z ~ 1, we use inste ad the func- 
tion presented in iFirmani fc Avila-Reesel (|2010h . who mod- 
ified them in order to obtain a continuous function from 
z = to z — 4. 



Each galaxy is assigned a certain cold gas mass 
by drawing values from a Gaussian distribution with a 
mean and standard deviation according to the stellar-to- 
gas inass_i^uiiction M (M « , z) (fitted to observations) given 
by IStewart et al, (2009). The extrapolation of this func- 
tion to very low masses ca n produce unrealistic gas frac- 
tions. To avoid such cases, IStewart et ail (|2009h proposed 
to assign an upper limit to the gas fraction given by 
/i im (z) = M g /M h \ Mm=3xl0 s M(S!Z such that M 9 (z,M») 
fn m (z)M h (z, M»). The value of M* = 3 x 10 s M Q was prob- 
ably chosen because is the minimum mass for which a mea- 
su rement of the c old ga s fraction was inferred in the work 
of IStewart et~afl l|2009ft . Given the high degree of uncer- 
tainty in the values of / gas for low-mass galaxies, increasingly 
higher for higher z, we opt to simply set /gas.ii m ^ 100, which 
is the maximum observed value reported in IStewart et al.l 
l|2009h . 

As a halo grows, the stellar and cold gas masses of 
its central galaxy are given by the aforementioned relations 
(with the new stellar mass added to the disc component), 
except during a period of time after final coalescence in a 
galaxy-galaxy merger. In this case, there is an "instanta- 
neous" increase of the stellar mass of the remnant that is 
naturally higher than the value given by the mean M*(Mh) 
relation. In such cases, the halo-to-stellar mass relation is 
only restored in the future once the halo has grown suf- 
ficiently. For some cases, this moment has not happened 
yet by z = and thus, the final sample of mock galax- 
ies has slightly higher stellar masses than predicted by the 
M*(Mh, 2 = 0) relation. We have verified that the shape of 
the stellar mass function that we recover at z = is roughly 
in agreements with that of the observe d relation by compar - 
ing to the best Schechter fit reported in lPanter et al.l {2007]). 
Overall, the values are always within a factor of ~ 2 — 3 with 
a spread that is larger for higher masses (because in average, 
more massive galaxies have more recent major mergers). 

The redshift we choose to seed the galaxy population 
is Zsced = 3 since for higher redshifts, the stellar and cold 
gas mass assignment becomes highly uncertain. Galaxies are 
seeded as pure discs at z see d, thus, we are effectively neglect- 
ing their previous morphological evolution. This is partially 
justified observationally since the fraction of disc-dominated 
galaxi es gets larger with h igher redshift. In particular, at 
z ~ l. lPannella et all (|2009l ) report a transition stellar mass 
of lO n M0 below which disc-dominated galaxies dominate 
the gala xy population . Simi larly, from extrapolating the re- 
sults of lOesch etakl (|201Gv i to z ~ 1, we estimate that 
~ 90% of all galaxies with M, < 10 n M o are disc-dominated 
(B/T < 0.5), whereas for larger masses, the mix is approx- 
imately 75% and 25% for disc- and bulge-dominated galax- 
ies, respectively. For z = 3, the mix would be even more 
in-balanced towards mostly disc-dominated systems. 



A2 Channels of merger-driven bulge growth 

A2.1 Transfer of stars from the primary and secondary 
into the bulge. 

In its orbit through the primary, the secondary looses en- 
ergy and angular momentum sinking to the center due to 
dynamical friction. In the last stages of coalescence, the sur- 
viving mass of the secondary M2 (dark matter + baryons) 
collides with the central region of the primary, of mass Mi . 
During this process, the collisionless components of both 
systems are subjected to rapid changes of the gravitational 
potential that broaden their energy distributions leading to- 
wards an equilibrium state. This violent relaxation process 
drives the stars originally rotating in discs towards random 
orbits forming a spheroidal remnant. The stars in the pri- 
mary disc affected by the dynamical action of the secondary 
are expected to be those within a radius enclosing the mass 
~ M2 . The stars at larger radii in the disc are also perturbed 
but likely they are re-arranged into final configurations that 
are not far from the original ones (for instance, one of the 
effects of the merger is to vertically heat the galactic disc). It 
is then reasonable to expect that approximately a fraction: 

fP _ m *,disc rcl _ M 2 . 

/rclaxcd " < disc ~ ^° ff = ~Mi (A1) 

of stars in the primary plus all the stars in the secondary 
will pass through violent relaxation and be transferred into 
the primary bulge. Due to the highly non-linear nature of 
mergers, they are better studied with numerical simulations, 
like those carried out in H09a. They seem to confirm the ap- 
proximation described before, although a more detailed -but 
uncertain- accounting of the efficiency of violent relaxation 
as a function of radius suggests a slightly non-linear depen- 
dence on /x e g . The corr ection factor to Eq. (|A1[I suggested in 
iHopkins et all (|2009bh is 2(1 + fi^)' 1 , with a = 0.3 - 0.6. 
We adopt this correction and use a = 0.3 for our fiducial 
model, noting that the values of /f claxcd are thus smaller 
than those given by Eq. I All making violent relaxation less 
efficient in the primary disc and slightly decreasing the final 
values of B/T, particularly at low masses. 

A2.2 Starburst from the transport of cold gas into the 
centre of the primary. 

A fraction /burst of the cold gas in the galaxies looses angu- 
lar momentum because the gravitational interaction during 
final coalescence generates a non-axisymmetric response in 
the galactic discs that morphologically resembles a bar. The 
resulting stellar and gaseous bars are however out of phase 
because gas is collisional and stars are not. Because of this, 
the stellar bar torques the gas bar draining its angular mo- 
mentum. In this way, the cold gas is effectively removed from 
the original discs and transformed into stars, during a star- 
burst, in the bulge of the remnant. This process is efficient 
within a region internal to a critical radius r cr it and thus, 
the fraction /burst (assuming that the discs have exponential 
surface density profiles) is given by: 

/burst = = A [l - (1 + r clit /r d )e { - r «^] (A2) 

where r^ is an effective gaseous disc scaling radius, and A is 
a normalization factor. The ratio of the critical to effective 
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scale radius depends on the merger mass ratio and relative 
orientation and orbit of the progenitors, as well as their stel- 
lar and gaseous c ontent. A parameteriz ation of this ratio is 
given in Eq. A2 of lHopkins et all (|2009bh . We find that using 
Eq. (|A2|) with A ~ f/0.26 gives a reasonable match to the 
simulation results shown in Fig. 7 of H09a. We note how- 
ever that the precise value of the normalization is not very 
relevant for our results (unless A>1). This is because the 
cold gas mass that participates in the starburst is usually 
subdominant compared to the stellar mass that is violently 
relaxed in the primary and secondary. 

The total mass of the stellar disc and bulge of the rem- 
nant after a merger are then given by: 



m *,disc(l /relaxed) 



, bulge 



,bulgc 



+ m », disc/relaxed + 7Tl * 



rp p i rs s 

+ Jburst m cold + /burst m cold 



(A3) 



where the superscripts p and s refer to the primary and 
secondary, and all the quantities on the right are defined 
just before coalescence; it is then necessary to estimate the 
dark matter and baryonic mass of the galaxies at this time. 

For the primary, the effective dark matter mass is that 
of the most bound material just prior to coalescence, which 
is well approximated by replacing the real density distribu- 
tion by a NFW radial density profile and computing the 
mass interior to a radius r s — r v i r /c, where c is the h alo 
concentration (the c(Mh,z) relation of iGao et al.l (|2008l ) is 
used). At this stage, its baryon content is simple given by 
the M*(Mh, z) and M 9 (M», z) relations defined earlier. 

The effective dark matter mass of the secondary is com- 
puted from the properties its halo had when it was about 
to enter the virial radius of the host. Due to tidal strip- 
ing, most of its dark matter mass will be lost before final 
coalescence, but the most bound material should survive. 
Because of this process, the subhalo mass in the final stages 
of the merger is an unreliable tracer of the potential well 
that shaped the properties of the satellite galaxy. For this 
reason, the empirical M*(M^) and M g (M t ) relations only 
apply to the secondary before it enters the virial radius of 
the main halo, and not afterwards as it spirals inwards. Mod- 
eling the evolution of the stellar and gaseous components of 
the secondary is an important element of the semi- empirical 
model we have used in this paper, we discuss it below. 



Initially, the disc has an an exponential gas surface density 
with an effective scale length given by 2 times the effective 
radius of the disc, which we take to be given by the following 
empirically motivated relatiorQ: 



r c s 



5.28kpc (1 + z) 



( M, 



(A5) 



V10 10 M O- 

We then apply Eq. (|A4jl for azimuthally averaged values of 
E 9 in rings to obtain the mass of new stars formed in a time 
interval dt as a function of galactocentric radius R: 



dPh(R) = (1 - Tl)-KA a {R)T, t (R)dt 



(A6) 



where A a (R) is the area of the annulus at radius -R, and 
1Z = 0.35 is the fraction of the mass in stars that is instanta- 
neously returned to the cold gas. This formula is reasonable 
as long as dt is smaller than the overall time scale of SF. 

For our fiducial model, we also adopt a recipe for a SN 
feedba ck-driven outflow following the energy-driven wind 
case of lDutton fc van den Boschl l|2009l) . The amount of cold 
gas ejected in each annulus is then given by: 

26efb-Esn??sn 



dM cjact (R) 



V£ C (R) 



-dM,{R) 



(A7) 



where Esn = 10 51 ergs is the energy of one SN, t^sn = 
8.3 x 10 -3 is the number of SN produced per solar mass 
of stars, and V eBC (R) is the escape velocity at radius R given 
by the gravitational potential of baryons and dark matter. 
The kinetic energy of the wind is assumed to be a fraction 
tEFB = 0.16 of the SN energy. 

Using Eqs. (|A4IIA7j) we estimate the gas and stellar con- 
tent of the satellite at the time of coalescence (recall that 
the merging time-scale is described in section |2~2")) . 

Since our model is not aimed to follow the bulge growth 
of satellites, we cannot assign B/T values to them. This im- 
pacts the amount of cold gas in the secondary that partic- 
ipates in the central starburst during a merger. By testing 
extreme B/T values for the satellite, from pure discs to pure 
bulges, we find that, in most cases, the morphology of the 
satellite is not a significant factor in the final morphology of 
the remnant. This is mainly because the mass that partici- 
pates in the starburst is usually lower that the stellar mass 
that is vi olently relaxed in the merging galaxies (see also 
Fig. 14 of lHopkins eralll2010al ). 



A3 Evolution of satellite galaxies 

A first simplistic approach would be to assume that the 
amount of cold gas and stars in the satellite just before co- 
alescence is the same as in the moment the merger started 
(the results in the first row of Fig. [3] are obtained using this 
approach). A more physically-motivated treatment, used in 
our fiducial model, is to assume that there is no further gas 
accretion (due to environmental processes such as starva- 
tion, ram pressure and tidal striping) but the satellite suf- 
fers SF and feedback proce sses. For the form er, we assume 
the Kennicutt- Schmidt law ()Kennicuttl il998l for gas surface 
densities E 3 above a critical density E c = lOM0pc~ 2 : 

E^BMepc-Gyr-^^l^)" (A4) 



APPENDIX B: VARIATIONS OF THE 
SEMI-EMPIRICAL MODEL 

We study now, in an almost qualitative way, the variations 
to some of the relevant ingredients of the semi-empirical 
model of galaxy occupation, concentrating on the statistical 
changes to the morphological mix of galaxies at z — 0. 

Changes in M*(Mh,,z). A substantial change in the 
evolution of the M*(Mh, z) relation has consequences for the 
efficiency of bulge growth thro ugh mergers . We e xplored this 
using the relation reported in I Yang et alj (|201ll ) for case 12 
in their Table 3. At z = 0, the difference betw een this for- 
mula and that o f F lrmam fc Avila-Ree"se] (|2010l ) is not signif- 
icant, except at the low-mass end where the former predicts 

5 Hopkins, P., p rivate communication, based on the results of 
IShen et alj J2003h . 
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lower values of M* for a fixed halo mass. At high redshift 
however, both formulae di f fer sig nificantly in the low-mass 
end, with the lYang et alj (|2011 ) formula predicting lower 
stellar masses for a given halo mass. This means that sec- 
ondaries will be less destructive for central galaxies. The 
net effect is an overall reduction in the values of B/T, ex- 
cept for the most massive galaxies, whose morphologies are 
typically shaped by massive secondaries. Thus, by using a 
M*(Mh, z) relation with systematically lower stellar masses 
in the low-mass end, the morphological mix changes towards 
a more extreme dominion of bulgeless galaxies in this mass 
regime. For example, for the case analysed, we find that at 
M* = 1O 1O M , ~ 85% of galaxies have B/T sC O.f , which is 
~ 1.7 times more than in our fiducial case, whereas there is 
essentially the same morphological mix for M, > 10 n M Q . 

Gas fraction limit. T he f sas (M*,z) empirical relation 
from IStewart et al. (2009) that we have used to assign a 
gas fraction to our mock galaxies can b e trusted up to z = 2 
down to M» ~ 3 x 10 8 M Q (see Fig. 2 of lStewart et al.ll2009h . 
For higher redshifts and lower masses, there is essentially 
no data to constrain this relation. In our fiducial model, a 
higher value of / gas at the start of the merger implies that 
the secondary will have more material available for SF as 
it spirals inwards towards the primary, resulting in a more 
destructive merger. We investigate the impact of a limiting 
value of /gas different to the one used in the fiducial model 
(/gas.iim = 100). This alternative value is given by the cosmic 
baryon fraction: M gas + M, ^ fbM^ with ft, = 0.17. In this 
case, all of our mock galaxies end up being more massive 
at 2 = 0, the change being significant for the least massive 
galaxies and negligible for the most massive ones. The values 
of B/T are systematically higher, specially in the low-mass 
end, which produces a morphological mix that is deficient 
in bulgeless galaxies. For instance, at 1O 9 M0, the models 
predict that only 20% of the galaxies have B/T < 0.1, while 
for M* ^ 10 M@, there is essentially no change compared 
to our fiducial model. 

SF and feedback recipes for the secondary. The satel- 
lite evolution that we have chosen in our fiducial model is 
based on the assumption that the Kennicutt-Schmidt law 
is valid at all redshifts.lt is possible however, that this law 
changes with z, being steeper and with lower normalizat ion 
at higher z (e.g., see Fig. 3 of lGnedin fc Kravtsovll2~010l ). A 
lower SF rate at higher z will result in satellites with higher 
gas fraction at the moment of final coalescence and thus, in 
mergers that are less destructive. We find that within the 
range of values of the normalization, exponent and threshold 
of the Kennicutt-Schmidt law at z = 0, the values of B/T 
do not change significantly when varying these parameters: 
slightly higher B/T ratios are obtained for higher values 
of the normalization and the exponent, and lower values of 
the gas density threshold for SF. Only a very extreme de- 
viation from the Kennicutt-Schmidt law at high redshifts 
would have a significant impact on the morphological mix 
of central galaxies at z = 0. Another ingredient in the satel- 
lite evolution is that of feedback-driven outflows, which are 
more effective for low-mass galaxies. If we remove this mech- 
anism, the B /T values increase in the low- mass end because 
SF is more efficient in this case. At 1O 9 M0, a model with- 
out feedback predicts a fraction of bulgeless galaxies of 45%. 
At masses higher than 1O 11 M0, there are no changes in the 
morphological mix. 



Finally, we note that in our fiducial model we have as- 
sumed that once the satellite enters the virial radius of the 
host, all its hot gas corona, is immediately stripped. This is 
not realistic for the most massive satellites that will be able 
to retain some of its gas for some time after accretion. If the 
cooling time of this hot gas is lower than the merging time, 
then the gas will condense to the center of the satellite and 
contribute to the cold gas mass and also to SF. Although 
without a proper model, it is not possible to know the im- 
pact to our results, we speculate that satellites, specially the 
massive ones, will likely be more destructive resulting in a 
morphological mix with more bulge-dominated galaxies at 
high masses. 
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